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I discuss recent applications of QCD light-cone sum rules to various form factors 
of pseudoscalar mesons. In this approach both soft and hard contributions to 
the form factors are taken into account. Combining QCD calculation with the 
analyticity of the form factors, one enlarges the region of accessible momentum 
1 transfers. 

> 

, 1. Introduction 

t-H ■ 

■ At the beginning of this talk let me quote Misha Shifman: 

, " Unlike some models whose relation to Nature is still a big question mark, 

Quantum Chromodynamics will stay with us forever." 1 
0^ . In the strong coupling domain, the quark-gluon gauge dynamics of QCD 

manifests itself in a form of hadrons. A comprehensive analytic description 
of hadrons and their interactions in QCD remains an unsolved problem 
that may also stay with us forever. Today we have at our disposal only 



?-h ' approximate methods and effective theories of hadrons, with an impressive 

progress achieved in the numerical simulation of QCD on the lattice. 

An approximate analytical calculation of hadronic observables has been 
made possible with the advent of QCD sum rules. 2 Not only the original 
SVZ method is still being extensively used, but also its "offspring", the 
light-cone sum rules (LCSR). 3 In what follows, I overview applications of 
LCSR to hadronic form factors. I will also discuss a possibility to enlarge 
the region of accessible momentum transfers by employing the analyticity 
of the form factors. 
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2. Hadronic form factors 

The simplest hadronic form factors parameterize electroweak transitions 
between two ground-state pseudoscalar mesons. A well-known example is 
the pion form factor generated by the quark electromagnetic current: 

«P + <z)Ij7>(p)> = (2p + g)^(Q 2 ), (1) 

where = ^u(x)'jf 1 u(x) — ^d{x)^i fl d{x) and Q 2 = —q 2 . This form factor 
remains a popular test ground of QCD models, and is being used nowadays 
to probe the AdS/QCD approach (see, e.g., 4 ' 5 ). 

An important constraint on any model of the pion form factor is the 
Q 2 — > oo asymptotics 6 ' 7 : 



9Q 2 



1 



(2) 



determined by the perturbative gluon exchange between the quark and 
antiquark constituents of the pion. In this factorization formula, the long- 
distance dynamics below the separation scale fx ~ \/~Cp is encoded in the 
the pion distribution amplitude (DA) (p v (u,fi). The latter is defined as the 
twist-2 part of the quark-antiquark vacuum-pion matrix element expanded 
near the light-cone: 

(7r(p)|«(x)[x,0]7„7i5d(0)|0) x2= o - -ip„U [ due lup - x ^(u, fi) , (3) 

Jo 

where [x, 0] is the gauge factor and \i <~ l/^\x 2 \. The pion DA is a universal 
object, in a sense that it also enters factorization formulae for other pion 
form factors. The knowledge of <p n (u, /i) is however not sufficient for a form 
factor calculation at finite Q 2 . A major question remains to be addressed: 
how large are the soft contributions to the form factor (1) which start from 
0(l/Q 4 ). 

Switching from the electromagnetic to the weak flavour-changing cur- 
rent uj^b, one gets the B — > 7r transition form factor: 

(n(p)\u^b\B(p + q)) = (2p + qU+„(q 2 ) + .. , (4) 

where ellipses indicate the presence of the second form factor, due to non- 
conservation of the current. The c — ► d, s weak transitions generate similar 
D — > 7r, K form factors. The role of large scale Q 2 is taken now by the 
heavy quark mass ttiq (Q = b,c). In these heavy- light form factors there is 
again an interplay of hard and soft quark-gluon interactions. One can write 
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an asymptotic factorization formula, similar to (2), containing f n (u, fi) to- 
gether with the heavy-meson DA. Again, such an ansatz is phenomenologi- 
cally incomplete, because soft contributions have to be included. They play 
an even more important role in heavy-to-light form factors, since they are 
not power suppressed (in 1/toq) with respect to the hard factorizable part. 

Heavy-light form factors are not only interesting hadronic objects. Their 
knowledge is of primary importance for flavour physics. For example, the 
B — > 7r form factor is needed to extract the quark mixing parameter \V u b\ 
from data on scmilcptonic decay B — ► nlui, whereas D — ► it and D — > 
K form factors are used to determine \V c d\ and, respectively, \V CS \. The 
exploration of current and future experimental data on exclusive B and D 
decays demands accurate theoretical predictions for these form factors. 



3. QCD light-cone sum rules 

The method of LCSR was described in detail in many papers, (see, e.g. 8-12 
and the review 13 ), hence I only give a very short outline. The central object 
is the correlation function of two quark currents between the vacuum and 
on-shell pion state. To give a familiar example, the calculation of the B — > 7r 
form factors starts from the following expression: 



d 4 xe ^ ^(7r(p)|^{^^(x)7^(^),6(0)^m 67 5<0)}|0)-F( g 2 ,(p + (? ) 2 )p Al + ...,(5) 

where only the invariant amplitude multiplying p^ is relevant for the form 
factor fg^(q 2 ). At q 2 , (p + q) 2 <C m 2 , the 6-quark propagating in the cor- 
relation function has a large virtuality and the product of 6-quark fields is 
expanded near the light-cone x 2 ~ 0, leading to OPE for the correlation 
function, schematically: 

FoPE(q 2 Ap + qf) =]T f d 4 xe^C {ra) (x 2 ,m b )(7r(p)\u(x)r a d(0)\0) 

i-i J 



d 4 xe iqx / dvC^ '(x 2 ,v,m b )(ir(p)\u(x)G a0 (vx)r b d(O)\O) , (6) 



where gauge factors are not shown, T a ^ are various combinations of Dirac 
matrices. In the above, the short-distance coefficients Crp a ) and C(r b ) stem- 
ming from the 6-quark propagator are calculated perturbatively, whereas 
the vacuum-pion matrix elements are expressed via pion light-cone DA's. 
The term with r a = 7^75 yields the twist-2 DA defined in (3) and, in addi- 
tion, the twist-4 DA's. The terms with T a = i-f5, <J^-f5 yield twist-3 DA's. 
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Furthermore, the soft gluon emitted from the 6-quark propagator enters 
quark-antiquark-gluon operators, shown in the second line in (6). These 
operators are decomposed in pion three-particle DA's of twist-3 and 4. The 
dominant twist-2 and twist-3 contributions are calculated to 0(a s ), obeying 
collinear factorization 9-11 . Subdominant twist-4 and three-particle con- 
tributions are power suppressed, justifying the truncated twist /Fock-state 
expansion. It is a future task to assess also the twist-5 and 6 terms of this 
expansion. 

Matching OPE and hadronic dispersion relation for the correlation func- 
tion at |(p+g) 2 | ~ fni,x, where x is an intermediate scale ^> Aqcd, one uses 
quark-hadron duality for the hadronic states above B meson, obtaining: 



Substituting in (7) the amplitude Fope calculated from (6), one obtains 
an approximate analytical answer for the form factor in terms of pion DA's 
convoluted with calculable coefficients. Inputs in this relation include the b 
quark mass (in the MS scheme) , a s , and the set of low- twist universal pion 
DA's. The decay constant Jb is determined from two-point QCD (SVZ) 
sum rule. For other details, see, e.g. 11 . 

Importantly, in the resulting LCSR for the form factor fg^, the leading 
twist-2 and 3 terms start from a 9 , with NLO corrections being 0(a s ) sup- 
pressed. Thus, in LCSR both soft (nonfactorizable) and hard (factorizable) 
contributions are taken into account and the soft one dominates. Impor- 
tantly, higher twist contributions are suppressed by (A-QCD/rrib) and/or 
(Aqcd/x)- Since the calculation is done at finite mj, a transition from b 
to c quark in (5) is straightforward, yielding LCSR for D — > ir,K form 



The universality of the method goes even further. Forming a 
vacuum— >pion correlation function of the currents j^ 71 and ujp^d, one 
obtains LCSR for the pion form factor F 7T {Q 2 ) 16-18 , valid at Q 2 > Aq C£) . 
This sum rule predicts a substantial soft contribution at intermediate Q 2 
and, simultaneously, reproduces the QCD asymptotics (2). Furthermore, 
employing vacuums baryon correlation functions with a baryonic inter- 
polating current and DA's, one is able to calculate various baryon form 
factors, e.g., the nucleon form factors 19 . 

The method of LCSR has certain limitations. First of all, there is no 
"direct access" to the hadronic form factor. The analytical expression in 
a form of factorized OPE is obtained for the correlation function, whereas 



oo 








factors. 8 - 12 ' 14 - 15 
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the form factor enters the pole term of the dispersion relation. One isolates 
it from the contributions of higher states to the dispersion relation, esti- 
mated using quark-hadron duality. This approximation introduces a sort of 
"systematic" uncertainty of the method. To keep it under control, one fits 
the threshold parameter Sq in (7) by calculating the B-meson mass from 
the same LCSR. A typical accuracy of the form factors, calculated by this 
method is estimated at the level of ~ ±15%, by varying all input parame- 
ters and scales within their adopted interval. One of the main uncertainties 
is the shape of the twist-2 DA (/^(u) expressed via Gegenbauer moments. 

The region of momentum transfer q 2 accessible for LCSR is restricted: 
Ib-k^I 2 ) 1S calculated at q 2 <C (m B — m K ) 2 , practically at q 2 < 12 - 14 GeV 2 , 
f^(q 2 ) at q 2 ~ 0, and F„{Q 2 ) at Q 2 >1 GeV 2 . It is therefore an important 
task to access other regions of q 2 where more data are available. This is 
possible, provided one uses the analyticity of the form factors. 



4. Employing the analyticity 

Hadronic form factors are analytic functions of the momentum-transfer 
variable q 2 . A typical dispersion relation which follows from the analyticity 
of the B — > 7r form factor: 

oo 

2{m B « - q z ) 7r J s — q z 

takes into account the singularities located on the positive real q 2 axis, 
the lowest one being the ground-state B* pole (with the B*Bir coupling 
defined as in 8 ). Starting from the threshold at q 2 = (tub + m^) 2 , there 
are branch points and poles, generated by hadronic continuum states and 
excited resonances with B* quantum numbers. Note that, due to the QCD 
asymptotics f Br[ (q 2 — ► oo) <~ 1/q 2 (similar to (2)), there are no subtractions 
in (8). Importantly, this dispersion relation is valid at any q 2 . Hence, a 
practical way to enlarge the accessible g 2 -region is to match the LCSR 
result to (8) at q 2 <§C m 2 B and analytically continue the dispersion relation. 
This however can only be done if a model/ansatz is introduced for the 
integral over the spectral density of higher states in (8), e.g., an effective 
pole 20 . 

A less model-dependent approach employs conformal mapping (for ear- 
lier uses see 21 ). One maps the complex g 2 -plane where (8) is valid, onto 
the unit circle \z\ < 1 in the plane of the new variable: z(q 2 ,t ) = 



v /t + -q 2 + ^t + ~to 



, where t+ = (tub + m n ) , and to < t+ is a parameter. 
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There are many applications 22 of this approach to B — > it and other form 
factors, combined with pcrturbative QCD bounds obtained from the uni- 
tarity for the 2-point correlation function. In fact, the bounds are usually 
not restrictive, and, since the maximal values of \z\ corresponding to kine- 
matical boundaries of B — > it or D — > n, K transitions are rather small, a 
simple Taylor series near z = suffices to parameterize the form factor. The 
last version of this series parameterization 23 advocates a truncated power 
expansion: 

f ^ iq2) = i - J/m 2 X%(V,io))\ (9) 

q * B * k=Q 

(with certain constraints on a/t), where the B* -pole near the threshold is 
isolated. Three or four parameters are sufficient for a reasonably accurate 
parameterization. Employing (9), it is possible to go beyond the region 
where a LCSR calculation is valid. For that one has to fit the coefficients 
afe to the LCSR result for the form factor in the "trusted "region of q 2 . After 
that, transforming q 2 — > z, one continues (9) over z beyond the initial region 
and finally transforms the variable z back to q 2 . 

5. Recent results for heavy-light form factors 
5.1. B — > tv form factor and \V u b\ 

Let me first quote our update of the B — > tt form factor calculated from 
LCSR 11 with the result: /^.(O) — 0.26+o.o3j where a recent very accurate 
determination of the MS 6-quark mass 24 was used. In this calculation no 
attempt yet was done to use the analytical continuation. On the contrary, 
in order to diminish the theoretical error for the form factor at q 2 = 0, the 
calculated form factor shape at < q 2 < 12 GeV 2 (the estimated region of 
validity of LCSR) was fitted to the q 2 -distribution in B — > nlvi, measured 
by BABAR collaboration. 25 This fit allowed to tighten the constraints on 
the Gegenbauer moments of the twist-2 pion DA , yielding : 

Mu, M) = 6u(l - u)(l + a 2 {ii)Cl /2 {2u - 1) + a^Cl' 2 {2u - 1)) , (10) 

with a 2 (lGcV) = 0.16 ±0.01, a 4 (lGcV) = 0.04 ± 0.01 (neglecting a 6 ,...). 
Note that these intervals are quite narrow and are within broader "world 
averages" (see e.g., 26 ). Deviation of the DA from its asymptotic form is 
thus rather mild. Finally, in 11 the form factor at q 2 = was used to extract 
\V u b\ from the data on B — > ttW. This and other recent \ V u b\ determinations 
from B — ► irlvi are summarized in the following table: 
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[Rcf.] 


calculation 


input 


\V uh \ x 10 3 


27 


lattice 




3.38±0.36 


28 


lattice 




3.55±0.25±0.50 


10 


LCSR 




3.5 ±0.4 ±0.1 


29 




lattice 8 LCSR 


3.47 ±0.29 ±0.03 


11 


LCSR 




3.5 ±0.4 ±0.2 ±0.1 


23 




lattice® LCSR 


3.54 ±0.24 



Importantly, both LCSR and lattice QCD determinations are in a good 
agreement with \V U \,\ — (3. 5^o 14) x 10~ 3 inferred from the CKM unitarity 
triangle fits 30 . Let me also mention a LCSR calculation 31 of B — > K and 
B s — > K form factors (needed e.g., for models of rare exclusive B decays). 



5.2. D — ► 7r, K form factors and \ V c d\, \V CS \ 

Following the same calculational scheme as for B — > it 12 and B — > K 31 
transitions, using the same set of pion and kaon DA's and MS c-quark mass 
from 24 , we calculated 12 the D — > ir and D — > if form factors. In the latter, 
the SU(3)fi violation is taken into account in 0(m s ). Predictions for these 
form factors from lattice QCD and LCSR are presented in this table: 



Method [Rcf.] 






Lattice QCD 33 

34 

35 


0.57 ±0.06 ±0.02 
0.64 ±0.03 ±0.06 
0.74 ±0.06 ±0.04 


0.66 ±0.04 ±0.01 
0.73 ±0.03 ±0.07 
0.78 ±0.05 ±0.04 


LCSR 14 

15 
12 


0.65 ±0.11 
0.63 ±0.11 

o.67±g:S? 


0.75 ±0.12 

7^+ 0A1 
u - ,o -0.08 



Our calculation was then used to determine \ V c d\ and \ V CS \ with an improved 
accuracy from the recent CLEO data 32 on semileptonic D decays, (see 12 
for more details and numerical results). 

In addition, the series parameterization similar to (9) was used to 
enlarge the accessible q 2 region. The form factors were calculated at 
q 2 < (still within the region of validity of LCSR), and fitted to the z- 
paramctcrization. The result was then analytically continued, to cover the 
whole kinematical region < q 2 < (mo — "V(if)) 2 of D — > tt{K)Ivi decay. 
The predicted shape of f^iq 2 ) normalized at q 2 = is compared with the 
one measured by CLEO 32 in Fig. 1. 
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6. Pion form factor 

As already mentioned, the LCSR approach allows one to calculate any 
hadronic form factor, provided the momentum transfer is sufficiently below 
the hadronic threshold in the channel of the transition current. It is there- 
fore interesting to return to the pion e.m. form factor F 7r (Q 2 ) which is quite 
sensitive to the pion twist-2 DA and recalculate it from LCSR 17 ' 18 . In this 
sum rule the twist-2 term with 0{a s ) corrections and twist-4,6 terms are 
taken into account. With the pion DA given by (10) and taking the remain- 
ing input from 18 , 1 recalculated the pion e.m. form factor at Q 2 = 1.0 — 5.0 
GeV 2 . The result is presented in Fig. 2, for the central values of the input. 
Conservatively, a ±15% uncertainty still has to be added. The agreement 
with the recent accurate data at Q 2 < 2.45 GeV 2 , obtained at Jefferson 
Lab 36 is encouraging. 

Furthermore, it is easy to write down an analog of series parameteriza- 
tion (9) for F n (Q 2 ) , which is simply a polynomial in z in this case (the 
lowest p-meson pole is far above the threshold) and continue the form factor 
to small Q 2 . The result shown in Fig. 2 is consistent with the direct mea- 
surement of the pion form factor in this region 37 and correctly reproduces 
HmQ2^ F 7T (Q 2 ) = 1 (within the accuracy of the sum rule). 

I conclude that the pion DA (10) used in calculating the heavy-light 
form factors successfully passes the pion form factor test. A more detailed 
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I 1 1 1 1 1 

1 2 3 4 5 

Q A 2 [GeV A 2] 

Fig. 2. The pion e.m. form factor calculated from LCSR 17 ' 18 at Q 2 = 1 — 5 GeV 2 
(solid), fitted to series parameterization and analytically continued to smaller Q 2 
(dashed), compared with JLAB data' 16 (filled points) and NA7 data 37 (open triangles) 



study will be published elsewhere. 

7. Photon-pion transition form factor 

(a comment added after the talk) 

The pion twist-2 DA also determines the QCD asymptotics 7,38 of the 
7*7* — ► 7T° amplitude: 

F r m = &l j ^'j^ e o(a.). (ii) 



where both photon virtualities, Q 2 = —q 2 and |rC 2 | ^ Q 2 , are sufficiently 
large. The 7*7* — ► ir° amplitude is however only measured 44-46 when one of 
the photons is almost real (k 2 ~ 0), in which case it reduces to the photon- 
pion transition form factor F^(Q 2 ) = F~<" w (Q 2 ,0). Note that at k 2 = 
the factorization formula (11) is incomplete, because the real photon has a 
long-distance, hadronlike component. Hence, substituting a certain model 
of tp w (u,/j,) in (11) and comparing the result at k 2 = with experimental 
data on F 77r (Q 2 ), is not yet a conclusive test of QCD. Calculations which 
take into account long-distance component of the photon can be found in 
the literature, (see, e.g. 39 ). 
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> 

CD 
CD 




15 20 25 

Q A 2 [GeV A 2] 



40 



Fig. 3. Photon-pion form factor Q 2 F 7 ""(Q 2 ) calculated from a combination of LCSR 
and dispersion relation^ , including 0(ct s ) corrections^ 1 , compared with the CELLO^^ 
(open triangles), CLEO^ 5 (full points) and BABAR Jf(> data (full triangles) 



A method to calculate F n '' K (Q 2 ), combining LCSR with the analytical 
continuation was suggested in 40 . The 7*7* — ► ir° amplitude was calculated 
using light-cone OPE, with both Q 2 = —q 2 and \k 2 \ large, adding also the 
twist-4 term to (11). (A calculation of the next, twist-6 term remains a 
future task.) Equating this amplitude to the dispersion relation in k 2 : 

V2f p F pn (Q 2 ) , 1 [°°ImF^(Q 2 ,s) 



F^ E (Q 2 ,k 2 )^ 



m„ 



k 2 



+ 



1 r 



k 2 



(12) 



"p Js 

where p and u contributions are combined in one resonance term with 
rrip ~ m u , one calculates the p — > ir form factor using the technique of LCSR 
and employs quark-hadron duality for the integral above the threshold so- 
After that a smooth transition to k 2 — > in the dispersion relation is 
possible, yielding the form factor F n " K (Q 2 ). The 0(a s ) corrections in this 
approach were calculated in 41 , and some recent updates can be found in 

42,43 

Intrigued by the new BABAR data 46 , I returned to the calculation of 
F 1 " K (Q 2 ) from (12), using the twist-2 pion DA (10), the remaining input 
from 40 , and adding the 0(a s ) corrections from 41 . The result is shown 
in Fig. 3 (for the central values of the input). At large Q 2 , as expected, 
Q 2 F 17! (Q 2 ) tends to a constant a . In spite of a reasonable agreement with 
the old CELLO and CLEO data at low Q 2 and with the new BABAR data 
up to q 2 < 15 GcV 2 , this calculation (as well as the one in 43 ) does not 



a Note that the 0(a s In 2 Q 2 ) terms in the NLO part of the form factor 41 and the effect 
of their rcsummation still have to be assessed. 
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reproduce the tendency of the increasing Q 2 F 77r (Q 2 ) at larger Q 2 visible in 
the data, although the experimental errors in this region are rather large. 
An additional measurement, e.g., by Belle collaboration is important to 
confirm these very interesting data. 
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